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DECLARATION OF PAUL W. DENT 

I, Paul W. Dent, hereby declare as follows: 

1) Not later than 1 1 February 2001, 1 conceived an invention tentatively entitled 
"Non-Polluting Pilot Code Radiation." The invention was conceived in the United States while I 
was employed at Ericsson, Inc. I prepared an Invention Disclosure describing my invention and 
presented it to the Ericsson Patent Review Committee on or about 16 February 2001 . A copy of 
the Invention Disclosure Cover Form is attached as Exhibit 1. 

2) Between 16 February 2001 and 20 March 2001, Ericsson implemented their 
normal invention review process. After the invention review committee made the decision to file 
the application, the IP department at Ericsson selected outside patent counsel to handle the 
application. Subsequent to this selection, Ericsson mailed an order letter to the selected outside 
patent counsel on or about 20 March 2001 , with a request to prepare and file a patent 
application. A copy of the corresponding instruction cover letter to the outside patent counsel is 
attached as Exhibit 2. 

3) Upon receipt of the order letter, outside patent counsel opened a new matter, 
assigned a patent attorney/agent to the matter, and sent a confirmation letter to Ericsson on or 
about 23 March 2001 . A copy of the confirmation letter to Ericsson is attached as Exhibit 3. 
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4) After meeting with me, outside patent counsel drafted a patent application. A first 
draft of the patent application covering my invention was mailed to me for comments on or 
about 25 June 2001 . A copy of the corresponding cover letter is attached hereto as Exhibit 4. 

5) Between 25 June 2001 and 13 July 2001, I reviewed the application and 
provided comments to outside patent counsel. After revisions to the application based on my 
comments, and preparation of at least one more draft, a final draft of the patent application 
along with formal paperwork was mailed to me on or about 13 July 2001 . A copy of the 
corresponding cover letter is attached hereto as Exhibit 5. 

6) I signed the declaration and assignment on 18 July 2001 and mailed them to 
outside patent counsel. Upon receipt of the signed paperwork, outside patent counsel filed the 
application with the U.S.P.T.O. on or about 26 July 2001. The filed application received U.S. 
Application Serial No. 09/915,895. 

7) The invention claimed in currently pending U.S. Application Serial No. 
09/915,895, of which I am the sole inventor, was conceived on or about 1 1 February 2001 and 
pursued with reasonable diligence until the filing of the corresponding application on or about 26 
July 2001. 

I hereby declare that all statements made herein of my knowledge are true and that all 
statements made on information and belief are believed to be true, and further that these 
statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under §1001 of Title 18 of the United States Code 
and that such willful false statements may jeopardize the validity of the application or any patent 
issued thereon. ^ 
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Inventor: Paul W. Dent 

NON-POLLUTING PILOT CODE RADIATION 



NOTE: THIS DISCLOSURE INCLUDES ALL THE MATERIAL FROM MY DISCLOSURE ENTITLED 
"MOBILE PHONE DEVICES USING COHERENT MACRODIVERSITY " AND FROM MY DISCLOSURE 
" INTERFERENCE CANCELLING CDMA RECEIVING SYSTEM « BECAUSE THE THEORY 
IS CLOSELY RELATED. NEW TEXT DESCRIBING THE NEW ASPECT IS ADDED AT THE END. 
IT WOULD BE SENSIBLE FOR THE SAME ATTORNEY TO HANDLE BOTH . 

1. FIELD 



The invention applies to the transmission of pilot codes from cellular 
base stations using macrodiversity transmission, in such a way that the 
pilot code radiation does not cause electromagnetic pollution of the spectrum 
which could degrade communications . 

EXEMPLARY CLAIMS : 

1. A cellular communications system comprising a network of base stations 
for communicating with a plurality of mobile subscriber terminals, 
compr is ing : - 

-A loopback system for feeding back at least a portion of signals 
received at said mobile terminals to said network of base stations; 

-A downlink channel estimation system for processing said fed back 
signal portions and forming estimates of the coefficients of multipath 
propagation from each of said base stations to each of said mobile 
terminals with which the base station is communicating; 

-A transmit signal generation system for forming signals for 
transmission to said mobile terminals based on information to be 
transmitted to respective mobile terminals and said coefficients 
and at least one other signal such that said mobile terminals 
receive substantially no interference from said other signal when 
said coefficient estimates are accurate. 

2. The system of claim 1 in which said downlink channel estimation 
system detects the amount of said at least one other signal in 
said fed back signal portions and thereby determines and corrects 
inaccuracy in said coefficient estimates. 

3 . The system of claim 1 in which said downlink channel estimation 
system processes said fed back signal portions to determine the 
amount of wanted and unwanted information signals including said 

at least one other signal received at each of said mobile terminals 



and corrects said coefficients estimates in dependence thereon. 



2. SUMMARY 



A network of base stations for communicating with a plurality of 
mobile subscriber terminals over multipath propagation channels employs 
coherent macrodiversity techniques in which more than one base station 
transmits a signal modulated with information for the same mobile terminal. 
The signals transmitted by the base stations are formed from information 
symbol streams intended for respective mobiles and a set of downlink multipath 
channel propagation coefficients. When the number of independent transmitting 
base station antennas employed is larger than the number of mobiles to which 
information symbol streams are transmitted, an information deficiency in being 
able to determine all the multipath channel coefficients from every base station 
to every mobile terminal is solved by including a number of pilot or dummy 
symbol streams in the formation of the transmitted signals, the number of 
pilot symbol streams being equal to the difference between the number of 
independent base station antennas and the number of mobile terminals served. 
The pilot symbol streams are added to the base station transmissions in 
weighted amounts based on downlink channel estimates such that the mobile 
terminals receive substantially no interference from said pilot symbol streams 
or from symbol streams intended for other terminals when said channel estimates 
are accurate. The mobile terminals transmit at least a portion of their 
respectively received signal to the network of base stations. After separating 
the mobile signals in the network using interference cancellation techniques 
a channel estimation system detects the amount of each symbol stream in the 
fed back signal portions including the amount of the symbol stream intended 
for each mobile, the amount of unintended symbol stream interference received 
by each mobile and the amount of pilot symbol streams received by each mobile, 
thereby determining error in the channel coefficients and correcting the 
errors in dependence thereon. 



3, BACKGROUND 



The prior art of diversity communications comprises using more than one 
alternate communications means either simultaneously or by adaptive selection 
to improve communications quality or reliability. The alternate communications 
means can employ either different frequencies, different time slots or periods, 
different antennas or different polarizations or any combination of the above, 
being the known respectively as frequency diversity, time diversity, 
space diversity or polarization diversity. 

The more usual form of diversity is receive diversity, in which the alternate 
communications means comprises alternate reception means for receiving the 
same transmittter , When the alternate reception means comprises a plurality of 
antenna elements forming an array, the outputs of which are coherently combined, 
the technique is also known as receive beamforming. When the multiple antenna 
elements are widely separated so that a signal received at different elements 
may fade in an uncorrelated fashion, even to the extent of being received at 
either one or another element and not both at the same time, diversity combining 
the antenna output signals cannot be equated to forming a narrow reception beam, 
but nevertheless has a large advantage for mitigating the effects of fading, 
and is termed "space diversity reception "• . These forms of receive diversity are 
exemplified in the following U.S. patents, which are hereby incorporated by 
reference herein. 



6,081,566 Method and apparatus for interference rejection with different 

beams, polarizations, and phase references 
6,006,075 Method and apparatus for transmitting communication signals using 

transmission space diversity and frequency diversity 
5,991,282 Radio communication system with diversity reception on a time-slot 

by time -slot basis 
5,963,874 Radio station arranged for space-diversity and polarization 

diversity reception 
5,878,093 Interference rejection combining with frequency correction 
5,867,791 Up link macro diversity method and apparatus in a digital 

mobile radio communication system 
5,499,272 Diversity receiver for signals with multipath time dispersion 
5,481,572 Method of and apparatus for reducing the complexity of a diversity 

combining and sequence estimation receiver 

Receive diversity according to the above patents relies upon knowledge of 
the relative phase and amplitude of the same signal received at different 
antennas. This may be determined for example by correlating the antenna outputs 
with known symbol patterns embedded in the signal for that purpose, known as 
channel estimation. Diversity reception using channel estimation is described 
for example in the above '272 patent. How to use the channel estimates in a 
computationally efficient manner to decode symbols based on joint processing of 
all the antenna signals is the subject of the above '572 patent. 



Transmit diversity is the less usual form of diversity, in which the transmitter 
uses multiple antennas to communicate to receivers each having a single 
receiving antenna. When the multiple transmitting elements are in relatively 
close proximity and form an array, the technique is known as transmit 
beamforming. Transmit beamforming is described in the following U.S patents 
which are hereby incorporated by reference herein. 



6,088,593 Two-way paging system and apparatus 
5,940,742 Two-way paging system and apparatus 

5,909,460 Efficient apparatus for simultaneous modulation and 

digital beamforming for an antenna array 
5,848,060 Cellular/satellite communications system with improved 

frequency re-use 

5,812,947 Cellular/satellite communications systems with improved 
frequency re-use 

5,631,898 Cellular/satellite communications system with improved 
frequency re-use 

5,619,503 Cellular/satellite communications system with improved 
frequency re-use 

5,594,941 A cellular / satellite communications system with generation 

of a plurality of sets of intersecting antenna beams 
5,642,3 58 Multiple beamwidth phased array 

5,594,941 A cellular / satellite communications system with generation 
of a plurality of sets of intersecting antenna beams 

In the above transmit beamforming patents, beamforming relies on the assumption 
that the relative phases and amplitudes of the propagation path from each 
antenna element to a given receiver are known, predeterminable characteristics 
of the array. In the above '941 patent to current Applicant, It is also 
disclosed that mobile receivers may provide feedback information to the 
transmitter using a reverse communications channel to assist it in the 
predetermination of these array characteristics. This technique is called 
"Mobile-Assisted Beamforming" . 

When however the propagation paths from each antenna to each receiver comprise 



multiple propagation paths of different delay, and each has an amplitude and 
phase that varies with time, it is not possible to use the above methods of 
coherent beamf orming . This usually arises when the multiple transmit antennas 
are widely separated. In that case, the prior art describes the non-coherent 
method of transmit diversity variously called SIMULCAST, MULTICAST , or 
TRANSMIT MACRODIVERS ITY . These terms generally apply to using more than one 
transmitting base station on different sites to transmit the same information 
to a receiver. For example, the civil ground-to-air radio systems -transmit from 
the ground to an aircraft using the three nearest ground stations. The three 
stations employ frequencies relatively offset by more than the highest audio 
frequency so that their beat notes are inaudible* Police radio systems also 
transmit using multiple base stations, but generally synchronise the frequencies 
so that the summation at the receiver of signals received from, multiple 
transmitters simply looks like the multipath reception expected for UHF 
ground propagation in any case, but provides immunity against one or other 
transmitter being shadowed from the receiver by terrain obstruction. When 
analog speech modulation is used, care must be taken to ensure that the 
modulation is in time sync and not relatively delayed between the stations. 
On the other hand, if digital transmission is employed, it can be advantageous 
to arrange that different transmitters transmit the same information but with a 
deliberate delay of +/- one digital modulation symbol. An equalizer in the 
receiver is able to combine these relatively delayed signals which exhibit 
uncorrelated fading, thereby achieving diversity gain. This diversity gain is 
not obtained if the signals are not delayed by one or more symbol periods. 
Using an equalizer for diversity combining signals received from two or more 
transmitters is described in U.S. patent no 5,327,577 {Uddenfeldt et al, 
reissued 2 Feb 1999) which is hereby incorporated by reference herein. 

Such non-coherent simulcast or macrodiversi ty transmission is also exemplified 
in the following U.S. patents, which are hereby incorporated by reference 
herein : - 



6,104,933 Method and apparatus for control of base stations in 

macro diversity radio systems 
5,930,248 Radio communication system selectively using multicast with 

variable offset time 
5,883,888 Seamless soft handoff in a CDMA cellular communications system 
5,845,199 Simulcasting system with diversity reception 

5,724,666 Polarization diversity phased array cellular base station and 

associated methods 
5,812,935 Cellular system employing base station transmit diversity 

according to transmission quality level 

CDMA cellular systems most commonly employ macrodivers i ty , which is also 
known as soft handoff, in which a mobile station on the edge of two cells 
receives the same information transmitted from both cells during the transition 
period between receiving service from one cell and later receiving service from 
the other cell, as described in the above '577 and '888 patents. Non-coherent 
transmit macrodivers ity for CDMA is also described in U.S. patent no. 5,940,445 
to Kamin. Kamin transmits a signal and a one-symbol/chip delayed signal from two 
antennas. Kamin furthermore does not associated one transmit power amplifer with 
one antenna, but employs the matrix-PA technique first described by Welti in 
U.S. patent no. 3,917,998 . The matrix PA technique provides the advantage of 
pooling the power of all the transmitters so that one antenna may transmit more 
than the power of one transmitter if required, providing that the sum of the 
powers transmitted by all antennas does not exceed the sum of the available 
transmitter powers. The advantage of the matrix PA is however unrelated to the 
advantage of transmit macrodivers ity, and in any case can only practically be 
employed when all transmitters are on the same site. The current invention 
however may be employed even when the transmitters- and antennas are not on the 
same site. 



Non-coherent transmit macrodiversity is less effective in some important ways 
than coherent beamforming. Coherent beamforming enhances the received signal 
both by the summing of the transmitted powers of all cooperating transmit 
antenna elements and by the focussing of the transmitted beam towards the 
receiver, providing antenna gain in addition. Macrodiversity however provides 
only the former and not the latter. Coherent beamforming has the additional 
advantage of diminishing the interference transmitted in other directions. 
This allows the same frequency channel to be re-employed in different directions 
without interference, increasing communications capacity, which can be measured 
in units of voice channels per megahertz per square kilomter for example. 
In this method, the receivers for which the same .channel is re-employed must be 
spaced by more than the beamwidth of the transmission beam. This can be many 
kilometers when the array of transmitting antenna elements is on the same tower. 

With the current invention of coherent macrodiversity transmission however, 
using antennas on widely separated sites allows discrimination between 
cochannel users that are as close as a few meters or even centimeters to 
eachother, enabling for example three cellular phones in the same room to use 
the same cellular channel simultaneously. 

Non-coherent macrodiversity by contrast increases the geographical 
spread' of the interference from a given signal and reduces the reuse of the same 
channel. It nevertheless increases capacity over the capacity that would have 
been achieved without macrodiversity, as many receivers are located in the 
border regions between different transmitter service areas, when a uniform area 
distribution of receivers applies, and thus benefit from macrodiversity, 

The research publications listed below have shown that wireless systems ^ 
using multiple transmit antennas and multiple receive antennas can achieve 
high datarate capacity. 

[1] E. Telatar, Capacity of Multi-antenna Gaussian Channels, 
AT&T-Be 1 1 Labs Internal Tech. Memo., June 1995. 

[2] G. J. Foschini and M. J. Cans, On Limits of Wireless Communication 

in a Fading Environment When Using Multiple Antennas, Wireless Personal 
Communications, vol.6, (no. 3), pp. 311-35, Klu-wer Academic Publishers, 
March 1998. 

[3] G . J. Foschini, R. A. Valenzuela, Initial Estimation of Communication 

Efficiency of Indoor Wire-less Channel, Wireless Networks 3, pp. 141-154, 
1997 . 

The above references are hereby incorporated by reference herein. The background 
to the techniques described therein is similar to that of the current invention. 
An NxN matrix describing propagation from each of N transmit antennas to each 
of N receive antennas is essentially invertable to produce N separated channels 
on the same frequency that can each be used to convey a basic datarate, thereby 
achieving N times the datarate in total. However, in this art, the aim is to 
increase the datarate to a single reeceiver having the N receive antennas, and 
is achieved through the ability to jointly process the output of the N antennas 
at a single location. Typically, a receive matrix multiplication is performed 
to separate the N channels, which are then decoded using N, single-channel 
modems. The receive matrix multiplication can however only be performed if 
the N receive antenna signals are available at the single location. 

In the current invention however, the N receive antennas are not colocated but 
may be miles apart, and each is associated with a different receiver. 

The current invention therefore aims to achieve similar advantage without the. 
increase of complexity entailed by equipping each with N antennas and receiver 
chains. 



Capacity is increased beyond that achieved with non-coherent transmit 
macrodiversity if the geographical spread of interference from a given signal 
can be diminished. The current application discloses coherent macrodiversity, 
which can be employed when the transmitters have knowledge of the multipath 
propagation channel coefficients describing multiple, relatively delayed 
propagation paths to each of a number of receivers that receive different 
information using the same channel from the same transmitters, herein 
abbreviated to Channel State Information or CSI. The invention ^ to be described 
provides the advantages of coherent beamforming, normally obtainable only with 
cosited antenna arrays, as well as the advantages of space diversity for 
overcoming terrain obstruction or uncorrelated fading, normally obtained only 
with widely spaced or non-cosited antenna elements. In the inventive 
coherent macrodiversity, the geographical spread of interference is reduced 
because the signals from multiple transmitters enhance one another by a 
maximum amount only at the intended receiver antenna location and do not 
enhance one another to the same degree in other locations even only 
centimeters away. 



3, DESCRIPTION 



The propagation from a transmit antenna to a receiver may be described by the 
following equation: - 

R(i) = Co.S(i) + Cl.S(i-l) + C2.S(i-2) + C^S (i-L+1) ...(1) 

where R(i) is the complex value of the received signal sample 
(with carrier frequency removed) at time instant i; 

S(i) is the symbol transmitted at time i and Cj is the complex 

number describing the phase and amplitude of the propagation 

path with j-symbol periods of delay (relative to the shortest path, Co). 

The combination of a signal with itself delayed in steps of one symbol or sample 
period may also be described by the use of the Z-transform. A discrete-time 
signal S that is delayed by one sample is denoted by z~ L .S, where 
"z" is the time advance operator and its reciprocal is the time delay operator. 

Thus the polynomial a + b z" 1 + c tT 2 times S means 

aS(i) + bS(i-l) + cS(i-2) . 

Thus equation (1) can be written R = C(z _1 ).S 

where C is now the polynomial, in z" 1 , with coefficients Co....C Irl , 
R stands for the whole received signal sample stream and 
S stands for the whole transmitted symbol stream. 

Denote by C^ k the polynomial describing multipath propagation from 

antenna k to receiver j , with Rj denoting the signal received at 

receiver j and T k denoting the signal transmitted by transmit 

antenna k. Thus the entire picture of propagation from all antennas to all 

receivers is described by the matrix equation 

*j = [c jk j.T k „. ... 



-1 



where [C] is now a matrix, each of whose elements is a polynomial in z . 

If we desire each receiver j to receive only an intended symbol stream Sj 
then the transmitted signals needed to achieve that are given by 

T k = [CJ^.S. 

Thus the problem at hand concerns how to invert a matrix [C] of Z -polynomials , 
and whether such an entity exists and is numerically well- or ill-conditioned. 

Matrix theory describes the procedure for computing the inverse of a matrix as 
follows: First, transpose the matrix then replace every element by its cof actor 
to obtain a matrix known as the adjoint matrix. Then divide each element by the 
determinant of the original matrix to obtain the inverse matrix. 

The determinant of a matrix is given by sums of products of its elements and is 
computable in a predetermined, systematic fashion. For example, for the 3x3 
matrix 



Cll C12 C13 
C21 C22 C23 
C31 C32 C33 



the determinant is 



Cll <C22 .C3 3-C32 .C23) -C12 ( C2 1 . C3 3 -C3 1 . C23 ) +C13 (C21.C32-C31.C22) 

The cofactor of element Cll is (C22 . C33 -C32 . 23 ) , which is therefore the first 
element of the adjoint matrix. The first element of the inverse matrix is thus 

(C22.C33-C32.c23) 
Cll (C22 .C33-C32.C23J-C12 { C2 1 . C3 3 -C3 1 . C23 ) +C 13 (C21.C32-C31.C22) 



When the C's are z-polynomials , it may be deduced that the numerator is a 
z -polynomial of twice the order of the C's while the denominator is a polynomial 
of three times the order. The inverse matrix comprises elements that have both 
a numerator z-polynomial and a denominator therefore. The signals to be 
transmitted must be processed by multiplication with respective inverse matrix 
elements and summed to produce the transmitted signals. 

There is no problem in processing a signal by multiplication with a numerator 
polynomial, which just involves delaying the signal by multiples of a symbol 
period, weighting and adding, i.e. an FIR filter. The signal might conveniently 
be assembled into blocks, such as TDMA bursts, and a whole burst processed at 
once within digital signal processing memory before being output to the transmit 
modulator circuits. The processed burst is longer than the original by the 
length of the impulse response of the processing z-function, but can be 
prevented from overlapping adjacent bursts by allowing a guard time between 
bursts, as is common in TDMA systems such as the GSM cellular system. 
In the case of a continuous signal such as CDMA for example, the signal can be 
segmented into blocks, processed in the same way as described above for a TDMA 
burst, and then the processed segments, which are now longer and overlap, 
are linearly superposed so that the extended tails add to the neighboring 
segments. It is acceptable for the impulse response tails to overlap 
neighboring segments so long as the intended receiver positions are the same 
for the neighboring segments. When the intended receiver positions are 
different for neighboring segments, the tails of a previous segment overlapping 
a new segment would cause interference to be transmitted in the direction of 
the new receivers. In this case, adequate guard times should be used, or 
alternatively a means to allocate receivers to timeslots based on position 



can be used so that the spatial pattern of receivers served in adjacent 
timeslots is similar.. 

Processing a signal by a denominator z-polynomial is more problematic however. 
Such a z-f unction is an Infinite Impulse Response (IIR) filter. The impulse 
response tails off exponentially as a function of the successive powers of 
the poles of the z-f unction, i.e. the roots of the denominator polynomial. 
Thus if a root has a magnitude less than 1, the impulse response decays; 
if the root has a magnitude equal to one, the impulse response rings forever; 
and if the root has a magnitude greater than one, the impulse response grows. 
The latter however may be dealt with by processing the hypothesized signal block 
backwards in time to apply those roots that are greater than unity. Thus, 
the signal block is passed backwards through an IIR filter comprised of the 
reciprocal roots instead, which are less than unity and therefore give a 
decaying time-reversed impulse response. The roots that were already less than 
unity form an IIR filter which is applied in the forward time direction. 
It remains to determine what to do about roots that lie close to unity. 
Due to their slow decay, excessive latency in processing the signal exactly 
would be entailed. If the signal is not processed accurately, by truncating the 
IIR response for example, there is a risk of causing excessive interference 
overlap with adjacent blocks because the exact inverse Omatrix would not have 
been applied. 

If the exact inverse C-matrix is applied to the information signals, the 
receivers will not only receive ONLY their intended information, but will 
receive it with multipath propagation cancelled already at the transmitter, 
since the resulting channel from the transmitter to the receiver is C.CT 1 , 
which is unity. 

Omitting division by the troublesome denominator polynomial is equivalent to 
multiplying each signal by the denominator polynomial, so the resultant channel 
from transmitting network to each mobile would be equal to that polynomial, 
which in the above 3x3 example was three times the length of the individual 
multipath channel polynomials C- k . The receivers can have an 
equalizer to decode symbols in the presence of normal amounts of multipath 
propagation delay, but perhaps not three times that amount. Therefore omitting 
the entire denominator may produce a net channel that exceeds the receiving 
equalizer's capability. However, one third of the factors in the denominator 
can be omitted and this is equivalent to multiplying by a polynomial of 
l/3rd the order of the denominator, thereby creating an artificial multipath 
channel only of the same length as the original multipath channel. 

According to one aspect of this invention therefore, the l/3rd of the 
denominator factors which are omitted are chosen to be those corresponding to 
roots of magnitude closest to unity, which are the terms that cause the most 
slowly decaying impulse responses as well as peaks in the frequency spectrum of 
transmitted energy. Preferably, the roots having logmagnitude closest to zero 
may be selected as the roots of magnitude closest to unity. 

In summary of the above improved system of communication therefore, the 
following is proposed: 

A transmit processor is provided with NxNxL channel coefficients describing the 
NxN multipath channels of impulse response length L to each of N receivers from 
each of an equal number N of antennas available for transmission. 

Treating each set of channel coefficients as a z-polynomial of order L and as 
an element of an NxN square matrix of such polynomials, the transmit processor 
forms the NxN adjoint matrix whose elements are polynomials of order (N-l)L,^ ■ 
and a determinant polynomial of order ML. 

The determinant polynomial is factorized to determine the NT roots and the 



L roots having smallest absolute value of log-magnitude are discarded, 
leaving (N-l)L roots forming a reduced determinant polynomial of order 
(N-l)L . The signals for transmission are then processed through FIR filters 
of order (N-l)L corresponding to the adjoint matrix polynomials and through 
an IIR filter of order (N-l)L corresponding to the reduced denominator 
polynomial as follows: 

Information blocks intended to be transmitted to different 
receivers are assembled in the transmit processor. Each information block 
is passed in the reverse time direction through an IIR filter formed from 
the reciprocals of roots of the reduced determinant that have a magnitude 
greater than unity and also in the forward time direction through an IIR filter 
formed from the reduced determinant roots of magnitude less than unity to 
obtain IIR-processed blocks. 

The IIR processed blocks are then FIR processed by matrix multiplication with 
the adjoint matrix polynomials to obtain transmit blocks. The transmit signal 
blocks may be further filtered to restrict transmission bandwidth and upsampled 
and converted to continuous time signals. The continuous time signals are then 
I/Q modulated on a radio frequency carrier, amplified in respective transmit 
power amplifiers and transmitted using respective antennas. Since IIR, FIR and 
bandwidth-restriction filtering are all linear processes, they may be carried 
out in any order. 

The above-mentioned polynomial lengths apply to the case where the three 
transmitting antennas are on the same sites or nearby sites. 

The modulators, transmit amplifiers and antennas need not be on the same site 
but can receive their transmit blocks from a common transmit processor, by 
optical fiber landline for example. Different sites must however synchronise 
the frequencies of their radio frequency carrier generators to the same stable 
reference, using GPS for example. The relative phases must be at least as 
stable as the changes in multipath propagation channel coefficients, which 
means that their relative frequency accuracy should be better than the Doppler 
frequency caused by relative motion between the transmitters and the receivers. 
This ensures that the mechanism used to keep the transmit processor provided 
with up-to-date and accurate channel coefficients is also sufficient to 
compensate for any phase difference or phase drift between sites. 
When the three transmitting antennas are on different sites far apart, 
the polynomial orders may be somewhat higher, as will be further discussed 
below . 



A system according to the above description is shown in figure 1. 
Information symbol streams intended respectively for receivers RX1,RX2 and RX3 
are input to numerical baseband processor (4). The processor (4) generates 
complex number streams having a real or in-phase part I and an imaginary or 
quadrature part (Q) representative of the information symbols modulated on to a 
carrier. Using knowledge of the multipath propagation polynomials Cij(z), the 
baseband processor then forms weighted combinations of the complex number 
streams and themselves delayed. The complex combination sample streams are 
then converted to continuous time waveforms using DtoA convertors and 
upconverted to a desired radio frequency channel using respective quadrature 
{I/Q) modulators (3A,3B, 3C) . The upconverted signals are then amplified" to 
desired transmit power levels using respective Power Amplifiers (2A,2B,2C) 
and transmitted by respective antennas (1A,1B,1C). The transmitted signals 
propagate to receivers RX1,RX2,RX3 by way of the* multipath channels Ci j . 
Receivers RX1 , RX2 , RX3 thereby receive radio frequency signals modulated 
only with respectively intended information symbols as originally input 
to baseband processor (4). In figure 1, antennas 1A, IB, 1C may be on the same 
mast or tower on the same site or even on different towers at different sites. 
In the latter case, the connection from the outputs of baseband processor (4) 
to the inputs of modulators (3A,3B,3C) may be made by auxiliary . intra-netwbrk 



communications means, such as optical fibre links, wireline or microwave . links , 
always observing the necessary precautions to maintain the phase coherency 
between the signals across antennas. Various means to accomplish this coherency- 
requirement are disclosed in the above-incorporated '941 patent to applicant. 

Figure 2 shows more detail of baseband processor (4) . 

Signals for transmission to respective receivers are first IIR filtered by 
the same ^IIR filter, in this exemplary implementation, the IIR filter being the 
determinant z-polynomial of the propagation matrix reduced in order by omission 
of factors corresponding to roots close to the unit circle. The number ^ of roots 
omitted may be anything from none up to the receiver's equalizer capability. 
Each omitted root requires the receiver to deal with one extra symbol period's 
worth of multipath time dispersion. If different receivers have different 
equalizer capabilities, the IIR filters 10A, 10B and IOC could be correspondingly 
different; with different numbers of factors omitted. The IIR filters may also 
differ by factors intended to optimize communication efficiency separately for 
each receiver, as will be further discussed below. 

After IIR filtering, the array of FIR filters (20) processes the signals. Each 
signal is processed by a corresponding row of FIR filters (20) . The output 
signals from FIR filters (20) are summed down columns, indicated by the + sign 
at the junction of the line from one output to another. The summed outputs are 
then the signals to be modulated and transmitted from respective antennas. 

In the case where the receivers are the same and IIR filters 10A, 10B . . . 10C are 
the same, the net channels by which respective signals reach respective 
receivers are identical, despite different propagation coefficients. In effect, 
the transmit preprocessing compensates exactly for the different channels. 
It is well known in the trivial case of N=l (one antenna transmitting to ^ one 
receiver) that precompensation for the propagation channel at the transmitter is 
not optimum. Such precompensation would use more transmit power to transport 
frequency components lying in a bad frequency region of the channel, having 
high path attenuation, and less power to transport frequency components in a low 
attenuation region. An optimum prefiltering would do the reverse, transmitting 
most of the transmitter power where the propagation was good and not wasting 
power where the propagation was bad. The known optimum solution for the single 
transmitter, single receiver case is called the "waterfilling" solution. 

When received signal power varies with time instead of with frequency, i.e. 
so-called flat fading, the prior art of CDMA systems such as IS95 teaches that 
rapid power control should be used to attempt to keep the received signal levels 
more or less constant. However, the waterf illing solution for determining how 
transmitted energy should be distributed over the frequency spectrum may be 
applied to temporally varying channels showing that it is better to 
transmit less power when the channel is bad and more when the channel is good. 
Temporal waterfilling is the opposite of the power control technique propounded 
in prior art CDMA systems. 

When received signal energy is a two dimensional function of both frequency 
and time, two-dimensional waterfilling should ideally be employed to determine 
the energy to be transmitted in each frequency subband at each instant. 

The waterfilling solution has the characteristic that no power is transmitted 
in regions of the spectrum where the received signal is instantaneously below 
a threshold signal to noise ratio. Producing a spectral null over ^ a broad 
spectral region however implies high order bandstop filtering, which may be 
impractical or produce time dispersion in excess of a receiver _ equalizer ' s 
capability to compensate. The current invention avoids these disadvantages. 

The principles for optimum prefiltering according to the invention 
for the multiple antenna and multiple receiver case will now be derived. 



Instead of describing the propagation channels by z -polynomials , they could be 
described by frequency response functions CjJc(w). The frequency response 
function can be obtained' 'from the z -polynomial simply be substituting 
exp(jwt) for z or exp(-jwt*) for z' 1 . The frequency response functions 
can be evaluated at each spot frequency to obtain a propagation matrix of 
complex numbers at the spot frequency. The matrix of complex numbers can be 
inverted by computing the adjoint matrix and the determinant to form the 
baseband processing required to transmit that spot frequency component^ of each 
signal to respective receivers in an optimum manner. This can be repeated for 
each spot frequency component to obtain the baseband processing required for 
each spot frequency component and thereby a set of inverse frequency response 
functions which form a matrix for pre-processing the transmit signals. 

The spot frequency case has also been formulated in the above-incorporated '941 
patent to Applicant for the case where the number of antennas (M) exceeds the 
number of receivers (N) , i.e for the overdimensioned case. 

If Cjk{w) is an N x M matrix of complex numbers describing the phase and 
amplitude of the path from antenna k to receiver j at a spot frequency^, and 

Sj (w) is the frequency component w of the information symbol stream Sj intended 
for receiver j, then the transmit signals Tk(w) that should be transmitted 
from antenna k in order to ensure that each receiver receives only its intended 
signal component, while minimizing the total transmit power needed, 
are given by: 

Tk(w) - C # [C.C # ]" 1 Sj (w) 
By evaluating C # [C.C # ] _:L at each spot frequency 

'w' a matrix of transmit preprocessing frequency response functions is found. 
These can be converted to a set of impulse responses by Fourier Transforming _ 
them if desired, so that the preprocessing of the signals Sj can take place in 
the time domain with FIR filters. However, it is possible that, in attempting 
to invert [C.C # ] at certain spot frequencies w, it is found that its 
determinant is zero or near zero. This would signify one of two problem 
situations: Either that there is an unavoidable null in the frequency response 
of the channel from all antennas to at least one of the receivers, or that the 
propagation channels from each of the antennas to two or more receivers are 
identical, causing the matrix to be rank-deficient. As in the z -polynomial 
formulation, zeros in the determinant of a matrix may be avoided by omitting 
that factor, corresponding now to a zero root, from the determinant- More 
specifically, roots in the frequency domain that lie on or close to the 
the imaginary axis represent "high-Q" poles and may be removed by multiplying 
the frequency response functions [ C . C # ] _1 C # 

with a zero term in the numerator having the same root. While this is in 
practice close to optimum, as in the time-domain formulation, this would result 
in all receivers experiencing the same net propagation channel, which seems 
unexpected and therefore likely to be non-optimum. Instead therefore, the 
following procedure is proposed: 

When the spot-frequency inverse matrix has been computed (the adjoint matrix^ 
elements and the determinant can be kept separate to avoid attempting to divide 
by a zero determinant), the net spot-frequency power attenuation factor frorn^ 
the M-antenna diversity transmitter to each receiver may be computed by summing 
the squares of the magnitudes of the adjoint matrix elements down columns. 
The squared magnitude of the determinant is then divided by the column sum to 
obtain a net power attenuation factor from the transmitting system to the 
receiver. Repeating that for all spot frequencies gives the power attenuation 
frequency response function from the transmitting system to each receiver. 
Mow the KT=1 theory can be applied to determine how a signal for a given receiver 
should be prefiltered for transmission through the asscociated channel 



having that determined net power/ frequency response. 

Prior art theory for N=l gives the optimum transmit power spectral shaping P(w) 
for transmitting information through a channel with a given power spectral 
shape H(w) as "the waterfilling solution". The waterfilling solution 
considers pouring water onto an object of shape 1/H(w) (with walls at the band 
edges} until the total amount of "water" contained is equal to the total power 
available. The function P(w) then equals l-Ho/H(w) where Ho is the recipro-cal of 
the level of the water. No power is transmitted in regions of the band still not 
submerged, i.e. where H(w) < Ho. The solution thus depends on the total power 
available relative to the receiver noise and the mean path attenuation, i.e. on 
the mean SNR that can be created at the receiver. Given a desired datarate and 
a channel bandwidth, the mean SNR needed to sustain that datarate can be 
determined, and the water-pouring solution for that SNR used. 

The water-pouring solution results in zero power transmission in regions of the 
spectrum where the SNR is less than a threshold. It is complicated to realise 
a transmit prefilter with abrupt cut-off at the boundaries of such regions, 
and these prefilters moreover have long impulse responses, requiring complex 
equalizers at the receiver. A compromise filter that centers single zeros 
of the form l+az" 1 on the null-regions can be used instead, which is 
nearly equivalent to the above-proposed time-domain solution of deleting 
determinant polynomial roots close to the unit circle. The difference is that, 
in the just-described approach, the position M a" of the zeros may be adjusted 
slightly away from the true determinant zeros in order to flatten the net 
frequency response from the diversity transmitting system to each receiver 
separately, taking account therefore of the effect of each column of adjoint 
matrix elements, thereby satisying the intuitive expectation that the net 
channels should not turn out to be identical for each receiver. 

The frequency-domain solution and the time-domain solution are in fact identical 
as the matrix 

C(w) # [C(w) .C(w) # ]- 1 

is just the matrix 

C(z" 1 ) # [C{z~ 1 ) .C(z" 1 ) # ]" 1 

with z~ l set equal to exp(-jwt), where the symbol # for matrices of 
z-polynomials means conjugate transpose and end-to-end (time-) reversal of the 
polynomials. An example of the correct way to form the time-reversed conjugate 
transpose of a matrix of z -polynomials is given in figure 11. 

Therefore the procedure can simply comprise determining the matrix C of 
z-polynomials as before, and then forming the matrix of double-length (2L) 
z-polynomials : 

Cfz" 1 ) .C(z- X ) # 

where # means time-reversed, conjugate transpose as defined in figure 13. 

The above matrix is inverted as before by computing the adjoint matrix of 
z-polynomials of length 2L(N-1) and the determinant z-polynomial now of length 
2LN, and then multiplying the adjoint matrix by C(z -1 ) # to obtain a matrix of 
length <2N-1)L polynomials. The attenuation/ frequency curve from the transmitter 
to each receiver can now be determined by letting z = exp(jwt) and summing 
down matrix columns the squares of the moduli of the elements, and repeating for 
each spot frequency. The square of the determinant modulus is then divided by 
the row sums. For each receiver, the resulting attenuat ion/ frequency curve may 
be examined to locate the highest L attenuation peaks. FIR zeros of the form 
-i 



1 + az" 1 are then placed on the attenuation peaks to either flatten 
the frequency response, which is close to optimum, or to over-flatten the 
response by turning the -peaks into troughs, which is even closer to the 
optimum water-pouring solution. An approximate rule-of -thumb for over-flattening 
would be to choose a zero to "turn a peak originally X decibels above the mean 
to a trough X decibels below the mean. The use of no more than L zeros ensures 
that the receivers do not need to equalize more than L symbol periods of 
net time dispersion. The result of applying the zeros to the signal will not 
however flatten the transmitter-receiver frequency response, as that is flat 
already when the zeros are not used; rather, it is the curve of transmit power 
versus frequency which is flattened across the spectrum, resulting ^ in notches 
in the frequency response from the transmitting system to the receiver. 
These notches appear to the receiver to have been caused by multipath 
propagation with a channel equal to the product of the applied "flattening" 
zeros . 

There is essentially no difference in the transmitter structure for the above 
solution as compared to figure 2, except for the choice of the FIR matrix 
polynomials and the IIR determinant polynomials, and prefilters 10A, 10B, IOC now 
also have an FIR element composed of the above-mentioned flattening zeros, 
if they do not exactly ' annihilate factors of the determinant. 
Note that if a determinant polynomial has a pole exactly on the unit circle, 
giving an infinite attenuation peak, then the pole would be annihilated by a 
flattening zero also exactly on the unit circle, i.e. |a|=l. Likewise poles 
close to the unit circle giving high attenuation peaks would be flattened by 
placing a zero nearby. If the zero was placed exactly on the pole, thereby 
annihilating it, the solution is identical to deleting the L poles of the 
determinant closest to the unit circle. 

Figure 3 is a flowchart for a control processor located in the network for 
determining the FIR and IIR transmit preprocessor coefficients based on channel 
state information (CSI). At step 100, the CSI information is updated to reflect 
latest estimates of the donwlink channel z -polynomials Cij . 

At step 101, a test is performed to determine whether the transmitter is trying 
to transmit' to a greater, lesser or equal number of receivers N compared to the 
number of antennas M at its disposal. If N>M for this implementation the 
requirement cannot be satisfied and an error must have occurred. Error recovery 
procedures at step 102c can for example comprise selecting to serve only the 
first M receivers and not the remaining N-M. If N^M, the matrix C is square 
and it is effectively inverted by computing at step 102a the adjoint matrix of 
polynomials of order (N-l)L and the determinant polynomial of order NL . 
Then at step 103 the determinant polynomial is factorized, i.e. its roots are 
found and sorted in order of magnitude. At step 104, the L roots and factors 
which are closest in magnitude to unity are discarded. The number L of roots 
and factors discarded may be different if different receivers have different 
equalizer capabilities, which the transmitter would be informed about in 
advance . 

If the number of receivers N is less than the number. of transmitting 
antennas or sites M, then the extra degrees of freedom are used to optimize 
performance at step 102b, which computes the matrix C.C # of 

polynomials of length 2L. Then the adjoint matrix of this matrix is computed as 
is its determinant. The adjoint matrix polynomials are of length 2L(N-1) and 
the determinant polynomial is of length 2NL . Other '.than that however, the 
procedure continues to the same step 103 as for N=M. After deleting the L 
roots closest to the unit circle, the determinant polynomial is now of length 
(2N-DL. Steps 105 and 106 are then performed with 'the resulting adjoint 
matrix and reduced determinant polynomial by the baseband processor of figure 2. 

If it is desired to perform optimum transmit spectral shaping for each receiver 
separately, the flow chart of figure 4 may be used. 



In figure 4, steps 103-105 of figure 3 have been repaced by steps 103a-105a 
of figure 4. At step 103a, the best achievable net attenuation versus frequency 
response from the transmitting system to each receiver is computed. It is then 
desired to warp the transmit spectrum for each receiver to direct more power to 
regions of the spectrum with more favorable attenuation and would have used 
less transmitter power and less to regions of the spectrum that suffer greater 
attenuation and would have used more transmitter power. The water-pouring 
algorithm is. the optimum and could be used at step 104a for this if not 
considered too complicated. Otherwise the simpler approach of adding flattening 
or overf lattening zeros can be used. In the case of high-Q determinant poles, 
i.e. roots very close to the unit circle, a flattening zero may be placed 
exactly over the pole to annihilate it. Instead of adding a zero then, a pole is 
annihilated from the determinant instead- On the other hand if one of the L 
poles closest to the unit circle is a low-Q pole, the attenuation frequency 
response may not show a peak exactly on the pole frequency but will be displaced 
due to the influence of the adjoint matrix FIR polynomials. In that case a zero 
is centered on the displaced peak and does not annihilate the nearby pole. 
Steps 105a and 106 are then performed by the baseband processor of figure 2. 



An example will now be given for the case of three transmitting antennas on 
different sites far apart. Figure 5 shows three base stations that 
for simplicity have been sited at the corners of an equilateral triangle with 
10km sides. The three sites are jointly supporting communications with multiple 
receivers, of which three are located at ml,m2 and m3 as shown, and it is 
desired to use the same channel frequency simultaneously if possible for all 
three in order to increase capacity according to this invention. 
The base stations transmit an exemplary symbol rate of 13MHz/48, which is the 
symbol rate of the GSM/ EDGE cellular system. The spatial spread of such a 
symbol is given by dividing the speed of light by the symbol rate, giving 
1107 meters. The distances between the mobile receivers and the base stations 
expressed in terms of symbol wavelengths is shown in the table below: 



RX1 
RX2 
RX3 



TXl 
3 . 10636 
5 . 58839 
4 . 65383 



TX2 
7 .1187 
5 .53673 
4 .66078 



TX3 
6 . 10841 
4 .56929 
6 . 67428 



If the channel for each of these nine paths were of equal amplitude and phase 
and were line of sight with no rnultipath, the z-matrix relating received signals 
Rl,R2,R3 to transmitted signals T1,T2,T3 would thus be: 



z 3 


■llRl = 


( 1 Z~ 4 


Z" 3 ) 

) 


Tl 


z 4 


1 

■ 56 R2 = 


( Z' 1 Z" 1 

r 


1 ) 

) 


T2 


Z 4 ' 


1 

,66 R3 = 


( 1 1 


Z" 2 ) 


T3 



where the common fractional symbol delay parts of each row and a common 
integer power of Z have been transferred to the left into R, such that the 
powers of z remaining are simplified to the lowest integers, representing 
whole symbols of delay . 

Maintaining the simplifying assumption of equal phase and amplitude on all 
nine paths, the adjoint of this matrix is: 



■ (Z" 3 -l) (V 3 -z~ 6 ) • o 

(l-Z" 3 ) (Z _2 -Z" 3 ) (z" 3 -z- 2 ) 

0 (Z~ 4 -l ) (Z^-Z" 5 ) 

and the determinant polynomial is -1 + z~ 3 + z~ 4 ~ z~ 7 = - ( 1-Z" 4 ) ( l~Z' x ) 

The determinant has all seven roots on the unit circle at 

2=1 (two roots) 
z - -1 
z - j 
z --j 

z = exp(jl20°) 
z = exp(j240 0 } 

Each root represents a frequency at which infinite attenuation can arise 
between the transmitting system and the receiving system, so it is inefficient 
to attempt to convey energy at those frequencies to the receivers. To avoid 
this, all seven roots on the unit circle should be annihilated by zeros in 
the numerator, which is the same as deleting the roots of the determinant. 

This could cause the receiver equalizer to have to deal with a length 7 
channel, i.e. 7 symbols of time dispersion. All the adjoint matrix 
elements however share at least one root with the determinant that can be 
annihilated. Cancelling the factor -1+Z" 1 from both the adjoint 
matrix and the determinant polynomial leaves: 



(l + z^ + z" 2 ) -Z~ 3 ( l + z _1 +z~ 2 } 0 
-d + z^ + z" 2 ) -z~ 2 z~ 2 

0 ( l + z _1 + z~ 2 + z~ 3 ) -z" 1 ( l + z -1 + z^ 2 + z" 3 ) 

and the determinant is now ( 1-z" 4 ) ( l + z _1 + z~ 2 ) 

Omitting to divide by the 6th order reduced determinant means that the receivers 
will receive their signals modified by a 6th order FIR filter, and their 
equalizers must be able to deal with 7 symbol periods of time dispersion. 

The above simplified example gave rise to many determinant roots on the unit 
circle due to the assumption of equal amplitude paths. If the paths are not 
equal, but exhibit a single propagation path of amplitude and phase given by 
the complex numbers Ci j , then the propagation matrix becomes instead: 
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The Cij may be complex numbers that are changing due to Rayleigh fading. 

If the paths are in addition multiple paths, the Cij can be z -polynomials as 

in the earlier examples. The additional powers of Z~ l that are now 

attached to the Cij account for the extra propagation delay due to the 

substantial path differences that exist in the case of widely separated sites. 
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A(7) = 
A{8)- 
A{9) = 
A(IO)- 
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A(12) = 
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REAL 
-0 .01492 
0 .01484 
-0 ,02419 
0 .00308 
-0 .04147 
0.36312 
0 .44006 
0 .58783 
0 .22755 
0 .33422 
0 .22015 
-0.67517 
0 .15489 
-0 .12943 
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ROOT (2) 
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ROOT (4) 
ROOT (5) 
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ROOT (ID 
ROOT (12) 
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0.01770 
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0.11521 

0.19125 

0 . 88261 

0 .01047 

1.09935 

-0 .36053 
0 .61808 

-0.26555 

-0.13633 

• al found by the above pro-am to be 

above ]wamG «,U»E \ 

0.46157 
0.01053 



0 



-0 . 



0 
-1 



-0 . 
0. 



REAL 
62535 
,35251 
21019 
47505 
,57 159 
.07257 
.76889 
.34753 
. 10402 
34591 
22513 
,95558 
43781 



-0 
0 
0 

1 . 



0 . 
-0 . 

-0 . 



-1 
-0 
-0 



IMAG 
07880 
94711 
.26791 
68848 
79850 
84329 
60238 
64730 
95889 
.29541 
. 10904 
.38249 
.97210 



DELETE 



1.0772E 
0.56192 
0 .018T6 
31070 
62521 
,30831 
.38001 
.29327 
0 .12369 
0 . 68946 
0 .96488 



0 . 
0 . 
0 . 

0 , 
0 



DELETE 



DELETE 
DELETE 



t to unity were determine^ by ^fofcLOG 
ots of magnitude closest U ^/ coffip lex logarithm func 



rr„ » — ; M b e dividing by tte ^ ■* 

have the nine rr 4 denominate 



Knlne regency responses snown ^ matrix 
have the nine tr q denominator and the ^ 

,v, P reduced determinant ^ ^/to obtain the nine-filter t 

deleted determinant facto 



closest to the unit circle; this effective channel comprises four rays that 
are exactly symbol spaced and contain the maximum possible energy. 

The nine frequency responses of figure 8 can also be combined in threes by 
adding their power responses to determine how much power in total is being used 
to transmit to each mobile receiver, as shown in figure 8C . The integral of the 
power spectral curves yields the total power used for transmitting the intended 
signals to each mobile. These powers can be compared to the powers which would 
have been necessary to communicate the same total signal power one mobile at a 
time from the best base station, with the same channel state to obtain gain or 
loss for the macrodivers ity system. In the above case, the results were : - 

Mobile 1: +4.46dB more power needed (a loss) 
Mobile 2: +0.98dB more power needed (a loss) 
Mobile 3: +2.50dB more power needed (a loss) 

The net loss is less than appears, because it is the mobiles that are relatively 
nearer their bases that receive the significantly increased power while those 
further away receive smaller power increases. 

After carrying out many trials using random propagation channels each 
comprising the exemplary three, symbol -spaced, delayed multipath rays, 
a histogram was obtained of the total transmitter power used compared 
to the transmitter power which would have been used for each base station 
talking to one mobile station independently, ignoring mutual interference. 
Figure 9 shows the histograms of gain/loss for this 3x3 system, 

the different curves corresponding respectively to the deletion of 1,2,3,4 or 5 
poles from the determinant. 

Figure 9A shows the gain/ loss histogram for the case where the number of 
poles deleted was selected in each trial to provided maximum communications 
efficiency. The mean loss in the 3x3 case is about 1 . 5dB in this case, for 
which the benefit of a tripling in capacity is obtained. It will be shown below 
that overprovision of base station transmitters in relation to mobile receivers 
can turn this net power loss into a net power gain. 

The logmagnitude of a root of the determinant polynomial indicates the rate of 
decay of the impulse response due to the associated root, in Nepers per symbol 
period. It could be decided that it was unnecessary to delete roots representing 
greater than 3dB {0.35 Nepers) decay per symbol period, for example. 

If only three roots had logmagnitudes less than 0.3 5, then only those need be 
deleted. Therefore an implementation of the invention can comprise deleting only 
those roots with absolute logmagnitude less than a threshold, and thus not 
always the maximum number that the receiver equalizer can handle, with the 
intention of approximating the performance of figure 9A for deleting the 
optimum number of roots . 

The key to coherent macrodiversity as described above is knowledge of 
downlink CSI at the transmitting network. This may be rendered possible under 
the following circumstances: 

(i) Uplink and downlink use the same channel frequency alternately in 
quick succession, a so-called time-duplex or ping-pong system. 
Then the transmitter may assume that the downlink channels are the 
same as it measures on the uplink when it decodes the signals received 
back from the receivers . 

(ii) The receiver measures downlink channel-related information, encodes it 
and transmits it back to the transmitting network with small turnaround 
delay. For example, the UMPTS Wideband CDMA system (W-CDMA) has the 
ability to serve up to 200 voice users per frequency channel per cell, 
or a proportionally lower number of high bitrate users such as 



mobile web-browsers. Therefore, for mobile web-browsers desirous of 
receiving a high instantaneous datarate, it is acceptable to use the 
whole capacity of a voice channel or more on the uplink to feedback 
channel state information . 

(iii) In a mobile satellite communications system where the relative 
coupling from transmit antenna elements to receivers is almost static 

\ depending only on receiver position. 

(iv) In a wireless-in-the-local-loop system for transmitting internet or 
voice services wirelessly to the home, and the receive antenna is fixed. 

(v) In mobile systems where the mobile terminal is likely to be stationary 
when high bitrate services are invoked. 

The above systems all provide feedback of channel state information, but the 
case of fast-moving receivers is the most challenging as the CSI changes 
rapidly, and low-delay, high-rate feedback of CSI is required. 
Some solutions for rapid feedback of changing CSI will now be described. 

Let [C] denote the current CSI assumed by the transmitter, which is 
in error from the correct CSI [C] by an error matrix [E] so that 

[C] = [C] + [E] , or conversely [C] = [C] - [E] 

The transmitter transmits [C'] _1 PjSj 
where Pj is the effective net channel 

for signal Sj . Pj is the factor by which prefilters 10A, 10B or IOC of figure 2 
differ from the true determinant polynomial. 

Receiver ( i ) receives 

Ri. - CC] ik . [CT^-jPjSj 

where summation over the common index k is implied 

= [C'-E^.tCT^PjSj 

= PiSi - [E] ik . [C rVjPjSj, 

- . 1 if i = j else 0 . 

since [C'] ik . [CT 1 ^ - d 1 ^ 

Thus if receiver (i) correlates its received signal with known symbols embedded 
in the transmission Sj to receiver(j), the error polynomial term 
[E] ik . [C ] k jPj summed over index k will be obtained. 

If all mobiles do this for all j including their own and return the results 
to the transmitter, since the transmitter knows the Sj it transmitted, the 
prefilters Pj it used, and the assumed CSI C f ij, it can compute Eij and hence 
correct C'ij towards the correct or changing Cij , thereby tracking changes 
in CSI. Thus the receivers can return measurements of uncancelled interference 
from symbols intended for other receivers, and delayed versions of the same. 
From these interference correlations, the transmitter deduces how its CSI must 
have been in error, and corrects it. Specifically, receiver 1 reports the 
polynomials determined by correlation with shifts of respective known symbol 
patterns as follows : - 

Xll(z) - PI - sigma{E lk C _1 kl} Pl 
X12(z) "= - sigma{E lk C" I k2} P2 



XlN(z) = - sigma {E lk C- 1 kN} PN 

This is a set of N equatipns for the N unknown polynomials El 1 , E12 , E13 . . . . E1N 

Likewise, receiver 2 reports 

X21(z) = - sigma{E 2k C _1 kl} Pl 

X22(z) = P2 - sigma {E 2k C" L k2} P2 



X2N(z) = - sigma {E 2k C" 1 kN} PN 

and this is a set of N equations for the N unknown polynomials 
E21, E22, E23 . . .E2N, 

and receiver N reports 

XNl(z) = - sigma {E Nk C _1 kl} PI 

XN2(z) = - sigma {E Nk C _1 k2} P2 



XNN(z) = PN - sigma{E Nk C" 1 kN} PN 
which is a set of equations for ENl , EN2 .... ENN . 

The solution of each of such sets of equations for one row of [E] is 

[C] [P_ 1 ] .X (1) 

where [P^] is a diagonal matrix of the reciprocals of the 

pref liters. If the reported measurements X were exact, the X polynomials would 
contain P as a factor which would cancel. The remaining factors would give a 
solution for E that was entirely FIR, i.e. no denominator polynomials, as 
required. However, due to noise, the reported X polynomials will not have this 
exact property. A solution is find the pure FIR solution of order L for E that 
best matches the frequency responses given by equation (1) . For example, 
denominator roots from P can be paired with the closest numerator roots from 
C or X for annihilation until only numerator roots remain. These then yield 
the "best" pure FIR solution for E. 

Figure 14 shows another method for providing channel state feedback from 
mobiles to the network. A network comprising network processing (20 0) and 
antenna sites (201,202) forms signals for transmitting to the mobile stations 
using the two sites. Mobile 100 comprises a transmit /receive antenna 101 
coupled via a duplexer 102 to receive circuits 103 which filter, amplify and 
convert the composite received signal to signal samples, preferably in digital 
form, i.e. using an AtoD convertor. The signal samples from the receiver are 
then added in 105 to a pilot code and then fed to transmitter circuits 104. 
The transmitter circuits convert the signal samples to an continuous signal 
using a DtoA convertor for digital samples, and the continuous signal is 
upcoverted to a transmot frequency, amplified to transmit power level and 
transmitted through duplexor 102 connected to antenna 101. 

The network sites receive the transmitted loop-back signal from various 
mobile transmitters. The loopback signals from different mobiles may be 
separated by interference rejection combining of the signals from the different 
sites in processor 200. Also, processor 200 computes correlations between the 
received signals and the pilot code inserted by a mobile to determine the uplink 
channel. Correlations are also computed between the network received signal 
from a mobile and the signals the network transmitted from each of its sites 



in order to determine the total loopback channel, which is the product of, the 
downlink and uplink channels. The uplink channel is then divided out to 
reveal the downlink channel. If necessary, the network sites can also each add 
a different, low-level pilot code to theor transmissions which would be chosen 
to assist in this loopback channel determination. Using this method, the 
mobiles are relieved of the complexity of performing channel determination. 



The mathematics given above include the derivation of the optimum transmit 
signals when the number of transmitters exceeds the number of receivers . 
With such overdimensioning, more degrees of freedom exist to find a way of 
overcoming fading on any particular path, and/ or of cancelling interference 
at each receiver. This should translate to fewer pathological determinant 
roots close to the unit circle, and therefore to greater transmit power 
efficiency. To obtain the benefits of overdimensioning it is only necessary 
for the number of independently fading paths of the same delay to be greater 
than the number of receivers. For meeting this criterion, mean delay from an 
additional base station can always be made the same as the delay from existing 
base stations artificially, which occurs automatically when the above-described 
mathematical operations are performed. However, this criterion can also be met 
by providing, at one or more base stations, transmitters for both polarizations, 
for example Right Hand Circular and Left Hand Circular. Field measurements 
confirm that the propagation channel for a transmitted wave of one polarization 
reaches a receiver over a landmobile propagation path shows little correlation 
to the propagation channel of the other, orthogonally polar ized .wave . 
Thus if each base station site Is equipped to transmit' power on both 
polarizations, the base stations each count as two independent base stations 
as regards the mathematical formulation, even though the RHC and LHC 
transmitters are co-located. 

Figure 10 shows a gain/loss histogram computed from many trials of a 
3 -base/ 3 -mobile system wherein each base station is equipped to transmit 
both polarizations. The channel propagation matrix is then a 3x6 matrix 
of polynomials, and the transmit signal generator computes six waveforms to be 
transmitted from respective base stations and with respective polarizations, 
based on the three information signals to be transmitted to respective mobiles. 
The different curves again correspond to the deletion of 1,2,3,4 or 5 of 
the "worst" determinant poles, and now show a net gain in total power 
efficiency of 4-4. 5dB in the mean. That means the sum of the power transmitted 
at all three sites and in both polarizations can be significantly less than 
needed in prior art systems not using the invention at the same time as 
allowing the same channel to be reused three times at each location. Moreover, 
the effect of using both polarizations to thereby overdimension the number of 
transmitters relative to the number of receivers is to virtually eliminate 
determinant poles close to the unit circle, so no determinant poles need be 
deleted. If no determinant poles are deleted, then the system provides a 
constant net propagation channel for information signals from the network to 
the mobile receivers, which is free of fading and multipath distortion. 

Strategies for slating determinant poles for deletion other than selecting the 
poles nearest to the unit circle can be used. If for example. the adjoint matrix 
polynomials in the same column all have a factor equal to or close to a 
determinant factor, then it is unnecessary to delete that determinant factor as 
it cancels with the same factor in the numerator. Likewise, when the numerators 
contain a zero close to an offending pole, they reduce the undesirable effect 
of the offending pole. A systematic way to detect and exploit such occurences 
to select the poles to delete for best effect is desired. This systematic method 
is provided using Cauchy's Residue theorem. Cauchy's residue theorem provides a 
much faster method to determine the total power transmitted when a "white" 
signal input is filtered by one or more numerator (FIR) z-po Lynomials and a 
denominator (IIR) polynomial to obtain the transmit signals. 
The method is outlined in an Appendix, and may be 



extended by a mathematician using the same principles to cover all practical 
cases. In the Appendix, a version of Cauchy's residue theorem is developed for 
integrating the modulus -squared of a complex Numerator frequency function 
divided by a complex Denominator frequency function. In the case where the 
Numerator is a z -polynomial at least one order less than the Denominator 
z -polynomial, the result is a sum of terms, each corresponding to the 
contribution from a respective denominator pole. If the numerator is small at 
that pole, then the contribution (the residue at the pole) to the total ppwer 
integral will be low. Therefore the systematic method desired above is 
to select for deletion poles that have the largest real part to their 
residues, inlcuding the numerator term when' computing the residues. 
It can be argued that the selection process should first select the pole having 
the largest real part to its residue for deletion, and then re-evaluate the 
residues of the remaining poles having deleted the pole. This process should 
continue until the desired number of poles have been deleted. The remaining 
denominator .order must be at least one higher than the numerator order. 
If the process shall be continued to denominators of lower order, then an extra 
residue appears due extra powers of of z on the denominator. The extra residue 
can be determined by making a partial fraction expansion to extract a term 
having that power of z on the denominator, the other term now having a 
numerator of lower order than the denominator. The extra residue is then equal 
to the coefficient of the highest power of z in the numerator above the 
extracted term. The extra term is needed to determine the power after 
deleting a pole causing the numerator order to equal or exceed the denominator 
order . 

In the application to this invention, the numerator can comprise the sum of 
the moduli of all the adjoint matrix polynomials down a column, where the 
modulus of a z-polynomial means the product of itself with its time-reversed 
conjugate, as further elaborated in the Appendix. The residue at a pole 
therefore comprises a part from each adjoint matrix numerator polynomial; so 
only if all adjoint matrix polynomials in the same column are small at the 
denominator pole with that pole have a small residue, and thus avoid being 
slated for deletion. 



When practising the invention, it is desirable to group mobiles using the same 
channel according to their locations repative to a group of three neighboring 
base stations. Figure 6 illustrates the difference 

between desirable and less desirable groupings. The desirable groupings in 
figure 6 are obtained by forming group 1 as comprising the three mobiles which 
nearest to their respective base stations. This can be considered as producing 
a propagation loss matrix with the least loss along the diagonal. Then group 2 
comprises the three mobiles with the second least loss to a respective one of 
the three base stations; group 3 comprises the three mobiles with the third 
lowest loss to their respective base stations, and so on. The least desirable 
groupings would comprise grouping together the three stations nearest to a 
single base station, for example. The reason for preferring the former grouping 
method is that it obtains the highest possible wanted-to-unwanted signal ratio 
even without employing the invention, such that the additional signals 
transmitted to cancel unwanted cochannel interference when practising the 
invention are a minimum. With the least desirable grouping, three mobiles 
receive cochannel signals from the same, nearest base station, and the other, 
more distant base stations must transmit a lot of interference-cancelling 
signal power due to their greater distance. 



It" is always desirable to simplify mobile receivers such as mobile cellular 
phones due to their large numbers, and place complexity instead in the network 
stations, which are much less numerous. Thus a simplified method by which the 
receivers can feedback downlink channel information to the transmitting network 



would be useful. For example, the signal received at each receiver could be 
simply turned around and retransmitted with minimum delay back to the network, 
as shown already In figure 14 . 

If all receivers do so* on the same channel, the network must separate them by 
uplink beamforming/ interference cancellation, which implies knowledge of 
uplink CSI. Uplink CSI is also needed to divide out the effect of the uplink 
channel polynomials on the retransmitted signal so that it reflects only 
the effect of the downlink channel. In a CDMA system, the mobile stations can 
retransmit the signal received on the downlink with the addition of an 
uncorrelated pilot code sequence that the network stations can use to derive 
uplink CSI . In a non CDMA system that would not tolerate an overlapping 
pilot sequence, the feedback signal can instead by periodically interrupted at. 
known times to insert pilot symbols that the netwotk stations can use to derive 
uplink CSI. Thereby the onus for analysing what the receivers have received is 
placed back on the network. The network has the great advantage of knowing 
every symbol that was transmitted to every receiver and what prefilters were 
used for all the signals. The network can therefore perform correlations using 
the entire symbol sequence transmitted to each receiver, including data symbols 
and not just known pilot symbols. Thus a modification to figure 14 can comprise 
interrupting the loopback signal to insert pilot symbols, replacing the 
additive combination of pilot and loopback signals formed by adder 105. 
In general, any suitable combination of the loopback signals with 
mobile-specific pilot symbols or mobile-discriminating information can be used. 

Many variations of the above principle of "mirror reflection" of the received 
signals back to the network can be devised. For example in a CDMA system, the 
received signal can be despread using the codes of each receiver to obtain 
despread symbols, then the despread symbols can be respread using corresponding 
uplink codes and added. The multi-code uplink signal is then mirrored to ^ the 
network. Interference correlations (the X polynomials in the above notation) 
can also be digitally coded of course, and transmitted as a data stream 
protected by error correction coding. For high symbol rates giving long channel 
polynomials (large L) or for large N (e.g. greater than 3) the amount of digital 
information to be transmitted may exceed the uplink capacity available, 
presumed to be for example the capacity of one voice channel, or about 
4 to 12 kilobits per second. The information could be selectively reduced 
by including in the reports only the X polynomial or polynomials having the 
greatest coefficient magnitudes; only polynomial coefficients that had changed 
by more than a threshold amount from a predicted value, or othe rmeans of 
down-selecting. Reporting only the coefficient with the greatest magnitude 
will cause the network to correct the transmitted signals to reduce only ^ that 
largest interference component, which however if repeated successively will 
sequentially reduce the interference in order of strongest components first. 

Another implementation of the invention with reduced CSI reporting requirements 
comprises the mobile stations reporting only the quality of the signal they 
are instantaneously receiving from the multi-antenna transmitting system. 
The transmitters transmit signals using a fixed or systematically varying 
phasing. The network examines the signal strength reports from the receiving 
stations and selects to transmit data intended for the receiver that is 
instantaneously reporting the highest quality increase over its mean received 
quality. Due to channel variations and/or the above-mentioned systematic 
variations in macrodivers ity phasing, it would be expected that each receiver 
received the downlink capacity for an equal fraction of the time. 

Thus, according to this embodiment, a group of for example three network 
transmitters transmit respective modulated signals derived from the same 
information symbol stream. The modulated signals are derived from the 
information symbol stream by varing the relative delay of the information symbol, 
stream for at least two out of the three modulation signals with respect to the 



third, such that, taking account of different propagation delay to the receiver 
for the dominant multipath ray, the modulation is received in synchronism at 
some receiver. Alternatively, the phase of the modulated radio frequency 
signal can be altered separately for at least two of the three network stations 
such that, taking account downlink channel phase for the dominant ray, the 
dominant rays are received in phase at some receiver. Preferably both phase 
and delay are altered for two out of three transmissions so that for some 
receiver, the dominant rays are received in synchronism and in phase. 

The receivers all demodulate the received signal and determine the quality of 
the demodulated data, which can for example be the total signal power formed 
from the sum of the squares of the ray amplitudes given by downlink 
channel estimation. The receivers report the measured signal quality 
continuously, for example by mirroring a signal related to the measured quality, 
combined with pilot symbol sequences to provide uplink CSI to the network. 
The network monitors the mirrored signal quality reports and chooses the data 
that it is instantaneously transmitting to be the data intended for the receiver 
intantaneously reporting the best signal quality relative to its mean reported 
signal quality. The number of receivers reporting signal quality in such 
implementations is not related to the number of transmitting antennas. 

One example of a macrodiversity transmission using the above principles would 
comprise three network transmitters, each transmitting the same information 
symbol modulation with no relative delay, but with relative frequency offsets 
of 0, -20Hz and +20Hz. The frequency offsets assure a varying net channel even 
if the receivers are static. 

Yet another implementation comprises a network of exemplary three neighboring 

transmitters. An unrelated number of receivers receives the signals from the 

three transmitters and each receiver reports the transmitter it receives 

most strongly and the signal strength. The network selects to transmit data to 

the receiver receiving the strongest signal relative to the mean reported 

signal strength, using the best transmitter. The need for full, complex 

downlink CSI is avoided in these implementations because no attempt is made 

to perform coherent interference cancellation in order to transmit three signals 

to three receivers at once. The above-described methods are different than 

the known method of choosing the destination receiver based on its 

instantaneous reception quality. These known method perforemd selection among 

receivers already assigned to a single base station, as opposed to the above 

methods which dynamically reassign the mobiles to receive from different 

base stations or conversely, dynamically reassign the base station that shall 

transmit to each mobile. 

In dynamically "re-sorting" the assignment of one of three base stations to 
one of three mobiles, there is a previously unreported gain in communications 
efficiency which is computed and reported below. 

In the prior art, a mobile station indicates its preferred base station and is 
assigned to that base station for a session. If fading temporarily hinders 
transmission, that base station is preferably temporarily used to transmit data 
to an unfaded mobile. The prior art does not attempt to do the above while 
using the same channel three times in an area served almost alike by three 
base stations. If this were attempted, the mean wanted-to-unwanted signal 
levels at each mobile (C/I's) may be computed from the mean distance loss 
using a 4th power propagation law, and averaged over a period long compared to 
the Rayleigh fading. On the other hand, when practising the current invention, 
the base stations reassign themselves to the three mobile stations dynamically 
using the feedback or loopback information in such a way as to minimize the 
total power transmitted for maintaining a unit received signal level at each- 
mobile. The average C/I's are then computed for the invention, and the 
following comparison obtained: 



Mobile 1 Mobile 2 Mobile 3 ■ 

C/I without adaptive allocation 4.47dB 4,93dB 5 . 12-dB • . 

C/I with adaptive allocation 4.81dB 5.31dB 5 . 14dB 

Gain for adaptive allocation 0.34dB 0.38dB 0 . 02dB 



The above gains were obtained with the relative base and mobile positions shown 
in figure 1. The mobiles received data from their nearest bases in the above 
tests 61.3% of the time, and fading caused re-allocation 39.7% of the time- 
It is interesting to place the mobile stations in greater proximity to eachother 
where the C/I would be expected to be very low with conventional fixed 
allocation. With three mobiles nearly equidistant from the three base stations 
and within 100 meters of eachother, the following results were obtained: 

Mobile 1 Mobile 2 Mobile 3 
C/I without adaptive allocation -1.65dB -1.61dB -1.62dB 

C/I with adaptive allocation 2,26dB 2 . 89dB 2.08dB 

Gain for adaptive allocation 3.91dB 3 . 50dB 3.70dB 



With three almost colocated mobiles, the allocation of base stations to 
mobiles stations exhibited almost equal probability of each of the six 
possible assignment arrangements. The arrangement that occurred most frequently 
during adaptive assigment was used as the baseline for the non-adaptive case. 

Thus there is a considerable mitigation of interference between almost 
colocated mobile stations when practising dynamic base station reallocation 
according to the invention. The invention guarantees a mean C/I of around 2dB 
minimum as opposed to approximately -2dB without the invention. 

This translates to a more than doubling of achievable data rates together with 
use of the same channel three times over in a given area. This gain may be 
achieved by reducing the spreading rate of a CDMA system for example, to take 
advantage of the C/I improvement by providing a higher data rate. Alternatively, 
the amount of error correction coding used can be reduced when the invention 
is employed, also increasing the useful information rate while maintaining 
acceptable error rates. The above advantages are obtained when practising fast, 
adaptive allocation during a transmission as opposed to the prior art of "slow" 
adaptive allocations made prior to the start of transmissions. Fast adaptive 
allocation switches base-mobile allocation in response to Rayleigh fading 
patterns, while slow adaptive allocation switches base-mobile allocations only 
in response to mean path-loss changes . 

In another implementation differing from the prior art, when a base station 
selects not to transmit to a faded mobile station, the faded mobile station 
may receive data instead from a base station to which it is unfaded, the network 
dynamically determining the base station to transmit the data based on the 
feedback signals. Alternatively, the network determines a mobile station and a 
base station in the vicinity that may momentarily use the channel without 
interference, based on which base-mobile link currently is the least faded 
relative to an equalizing threshold. The equalizing threshold is chosen for 
each mobile station to provide approximately equal probability of each mobile 
station receiving data, and therefore equal mean data rate. 

■A degenerate case of the M-antenna, M-receiver solution is the use of two 
antennas, for example two colocated RHC and LHC transmitting antennas, to 
transmit a signal to a single receiver. When the two antennas are colocated and 
differ only by polarization, the channel coefficients for the multipath channels 




are substantially uncorrelated for the two polarizations; however, the 

delays of the fading, multipath rays are the same. That is, if the channel for 

RHC is expressible as 

C^Z) = c o mc + c 1 RHC z- pl + c/^z-^ 2 

then the LHC channel is expressible as 

c lhc (z) = Cq LHC + Ci LHC z -pi + c 2 LHC Z" p2 

The power of Z are the same for both polarizations, expressing equality of the 
ray delays, but the coefficients are different and fade in an uncorrelated 
fashion. 

The optimum way to transmit a wanted signal to the receiver is then given as 
before by 

C # [C.C # ] _1 

times the wanted signal waveform, where C is now the 2x1 matrix of 
polynomials 

{ C™ c ) 
C = { ) 
{ C LHC } 

ThUS [CC*]" 1 - 1/(C RHC <C #RHC + C LHC_ C #LHC) 

and there is no adjoint matrix. The implementation then depends only on the 
location of the poles of this denominator /determinant polynomial. It may be 
appreciated that, when PI and P2 are large integers, this determinant has a 
prodigious number of roots even though it has a small number of coefficients. 

The above channels, which have three significant multipath rays of relative 
delays 0, PI and P2 symbol periods, are typical of a CDMA channel. 

A known CDMA receiver would conventionally employ a form of multipath equalizer 
known as a RAKE receiver, which correlates the received signal with a locally 
generated despreading code shifted in time by 0, PI and P2 symbols. 
The correlation results are then weighted and added. The number of delays 
used and correlations combined is called the number of " RAKE taps V The RAKE 
taps in the present example are placed at delays of 0, PI and P2 symbols 
or chips. The RAKE receiver is not limited so much by the amounts Pl,P2 of 
delay, but in the total number of taps. Thus, if a three-tap RAKE receiver were 
provided for use without the inventive coherent macrodiversity transmit scheme, 
it would be desirable that no more taps should be needed with the invention. 
However, omitting the whole of the above determinant polynomial would result in 
a net transmission channel having seven delays and needing in principle a 
7 -tap RAKE receiver with tap delays of 

-P2, -Pi, -(P2-P1), 0, {P2-P1}, PI and P2 chip periods. 

If the already-described pole-deletion procedure is used to improve the 
conditioning of the denominator for implementation, only two roots may be 
deleted while preserving only a 3 -tap channel for the RAKE receiver to deal 
with. Moreover, the three tap channel would always comprise delays of exactly 
0 , 1 and 2 chips with this method. Thus the three RAKE, taps would be adjacent 
instead of at 0, PI and P2 symbols delay. However, the RAKE receiver is limited 
only in the number of taps available to combine, and not their spacings . 
Therefore any three -coef f icient numerator polynomial of the form 



AO + Al . Z Q1 + A2 . Z Q2 

could be used as an FIR prefilter, and would preferably share roots with the 
denominator that it is desired to cancel. Alterntively, a polynomial in which 
most of the energy was contained in three taps could be used. 

A denominator polynomial of the form C 1 (z><: 1 # {z) + C 2 (z}C 2 # (z) +■ 

has the property that half of its roots are the conjugate reciprocals of the 
other. half. Those roots having magnitude greater than unity { and which would 
therefore require time-reverse IIR processing) may be slated for deletion. 
The remaining poles form a causal IIR filter which may be implemented 
entirely by forward-time processing, and the deleted poles give rise to a net 
FIR channel. This is one way of factorizing the denominator polynomial into 
the form C equ ( z ) C egu # { z ) , which factorization 

is always possible by selecting one root to form C equ and its 
conjugate reciprocal to form C equ # . The choice of the 
larger or smaller root for C equ can be made in 2 N/2 

different ways, given N/2 pairs of roots, so the number of such factorizations 
that can be made is large. 

Trial computations showed that factorization into a purely causal and a purely 
anti-causal part gives a channel/ although not limited to the same number of 
taps as the original channels CI and C2 , nevertheless contains most of its 
energy in relatively few coefficients/ allowing a RAKE receiver with a small 
number of taps to be used. When fewer taps are used than the number of effective 
channel rays containing signal energy/ a loss of energy is entailed. 
Figure 12 plots the PDF of this energy loss for using 1 to 6 taps to receive 
the signal through the above-mentioned anti-causal factor. The RHC and LHC 
channels were both of the form 

Co + CI . z 5 + C2 . z 13 / 

and thus the determinant is of order 26. The causal and anti-causal factors 
are then both of order 13, but do not have non-zero z-powers of only 5 and 13. 
Nevertheless, 80% of the energy is most often contained in the tap delay 14, 
93% is contained in tap delay 14 and tap delay 9 and 98% is most often 
contained in tap delays 14, 9 and 1. Thus it appears that a RAKE receiver with 
no greater a number of taps than would have been used for the conventional 
channel will be more than adequate for use with the Invention. 
Figure 13 plots a histogram of the net gain in communications efficiency 
for this method/ when the number of RAKE taps is fixed at three, the same as 
the number of rays in the multipath channels. The gain of 2dB equates to 
a 60% increase in communications capacity in terms of number of users or 
datarate when practising the invention. 

Figure 15 shows another application of the invention to a CDMA system. 
In this case, only a single transmit antenna of a single polarization is 
available to transmit to N mobile receivers. 

The signal to each mobile is spread-spectrum coded using a unique spreading 
code, and the set of spreading codes may be mutually orthogonal, at least when 
time-aligned. When not time-aligned, such codes will show some degree 
of correlation or mutual interference. 

The conventional CDMA system comprises inputting symbol streams S L . . . . S N 
to spread spectrum encoder 3 0 0, where is the coded 
information stream intended for mobile (i). Stream is multiplied 
with spreading code C and transmitted to mobile (i) through downlink 



channel polynomial P.(z). 

Mobile (i) correlates the received signal with a first shift of code C- 
and obtains Po(i).S i where Po ( i ) is the first, coefficient of P^z) and 
likewise obtains Pl(i).S i for other code shifts. 
These are then RAKE combined into a single sample stream whose 
dominant component is then, hopefully, Q°ii-S i 

where Qo iL = |Po| 2 +|Pl| 2 ,.|p l _ x | 2 

and L is the length of the downlink channel to mobile (i) . 

However, the combined RAKE output signal will contain other terms. 
There will be some Intersymbol Interference from the previous and the 
next symbol of stream S^^ due to correlation of mobile (i) code 
shifted with itself. If the channel time dispersion is less than one 
information symbol in duration however, only one previous and one future 
symbol will contribute to the RAKE output. This can be accounted for by 
coefficients Q_ 1 and Q +1 of a 3 -coefficient z-polynomial 
Q ii (z) and the combined RAKE output can be written 

Q u (z) .S i [z) 

In addition, there are correlations between code(i) of mobile (i) and code(j) 
of mobile(j) shifted, which give terms j ( z ) . Sj { z ) . 

The order of may, as already seen in the above simplified 

example, be extended when signal (j) is received from a base station lying at 
a greater or lesser distance from mobile (i), but not when all signals are 
received from the same base station. Thus all Q — in this example 
have three coefficients. This matrix of z -polynomials described the coupling 
from other mobile's intended information symbols to a given mobile's RAKE 
output, amd is nmot precompensated at the transmitter in the prior art. 

Using the invention however, this matrix of z -polynomials may be inverted 
as before by forming the adjoint matrix and the determinant polynomial. 
If the determinant polynomial has no unpleasant roots, then using the 
inverted matrix as precombiner 301 cancels all interference at the 
receiver output of each mobile from both other mobile's symbols and from its 
own previous and future symbols. If unpleasant roots of the determinant 
polynomial have to be deleted, and the remaining, reduced determinant used 
as the IIR part of precombiner 301, then all interference from other mobiles' 
signal will still be cancelled, but there will be a net multipath channel 
created from the product of the omitted determinant factors through which each 
mobile will receive its intended symbol stream. This can be completely 
compensated at the mobile receivers by following the RAKE combiner output with 
a 3 -tap, Viterbi Maximum Likelihood Sequence Estimator or other known form of 
symbol-spaced equalizer running only at the information symbol rate and not 
at the spread spectrum chip rate. Thus figure 15 contemplates the addition 
of an information-rate matrix precombiner to combine information symbol streams 
sl,s2...sn to form information-rate sample streams S1,S2....S N 
which are input to the conventional spread spectrum transmitting device 300. 
The precombiner 301 has the same structure as already shown in figure 2. 
The polynomials used in precombiner 301 are computed using knowledge both of 
the downlink channels P t ( z ) and of the cross -correlations between and 

and the autocorrelations of the different spreading codes used in coder 3 00, 
and using the assumption that the mobile receivers employ conventional RAKE 
combining with weighting coefficients adapted to the chip-level channel 



.impulse responses. 

Coherent macrodiversi ty transmission as described above relies on a knowledge 
of the downlink multipath channel coefficients in order to form the signals 
transmitted by the base stations. Channel coefficients are generally estimated 
by correlating the received signal with the transmitted signal to determine how 
propagation modified the transmitted signal to obtain the received signal. 
If the mobile subscriber terminal performs such channel estimation, the 
base station transmission generally includes symbol patterns called 
"a sync-word", "pilot symbols" or "pilot codes" known in advance to the mobile 
terminal with which to perform correlations. However, the base station network 
knows all the symbols transmitted, and is therefore in a better position to 
perfrom correlations using all transmitted symbols as known symbols, providing 
the signal received by the mobile terminal is available to the network. Thus the 
optimum method to determine downlink channel estimates is by feeding back from 
the mobile terminals to the network a copy of the signals received. 
Such loopback is possible if the mobile terminals are principally desirous of 
-receiving information from the network rather than desirous of transmitting 
information to the network (e.g. in mobile web-browsing applications) . 
When the mobile terminal loops back received signals to the network, the network 
will receive the looped back signals propagated also through the uplink channel. 
Thus the looped back signals received at the network are equal to the signals 
transmitted by the network propagated through the product of the uplink and 
downlink channels. In order to divide-out the effect of the uplink channel, the 
mobile terminals should therefore preferably add a syncword, pilot symbols or 
pilot code known in advance to the network with the aid of which the network 
can estimate the uplink channels, as was already shown in figure 14. 
When however the number of independent base station antennas N used to form 
macrodiversity transmissions exceeds the number of mobile terminals M, there is 
an informational deficiency in the network being able to determine all MN 
downlink channels from the loopback signals. This is most easily seen in the 
simple case of transmitting one information stream using two antennas 
(e.g. Right and Left Hand Polarization) to a single mobile terminal, as shown 
in figure 17 . 

Figure 17 illustrates a chip, symbol or signal sample stream X(z) being 
formed using a transmit preprocessor into a RHC transmit signal Pl(z).X(z) 
and an LHC transmit signal P2(z).X(z) where Pi and P2 are prefilters 
determined by and therefore known to the network. These transmit signals 
propagate to the mobile terminal through an RHC channel CR ( z ) and an LHC channel 
CL{z) respectively, yielding a signal Y ( z ) = ( PI ( z ) . CR ( z ) +P2 ( z ) . CL ( z } ) . X ( z ) 

By correlating the loopback signal with X(z) and correcting for the uplink 
channel the composite channel Pi ( z ) . CR ( z ) +P2 ( z ) . CL { z ) is obtained. However, 
this cannot be uniquely partitioned to determine CR(z) and CL(z) given Pl(z) 
and P2(z), as witnessed by the fact that the expression is identical to 

Pl(z) (CR(z)+aP2 (z) ) + P2(z) (CL(z) -aPl(z) ) 

since the terms Pl(z).aP2(z) - P2(z).aPl(z) cancel. Therefore it is not possible 
to determine CR within an arbitrary multiple 'a' of P2(z) nor CL within an 
arbitrary multiple 'a' of Pl(z). To solve this according to the invention, 
another information stream or pilot code must be transmitted jointly by the 
same base stations, but in such a way that the mobile terminals receive 
nominally no interference from. the pilot code. In prior art macrodiversity 
methods typically used in CDMA systems, each base station transmits a unique 
pilot code without attempting to prevent its reception at the mobile terminals; 
on the contrary, it is intended to be received by the mobiles and used to 
determine the selection of base stations to be used for macrodiversity. 
In the prior art, the total interfering energy of the pilot codes of all base 
stations can cause significant capacity loss, and is referred to as 
"pilot pollution". This degradation is avoided when practising the invention. 



The solution for the dual-polarization, single mobile terminal example of 
figure 17 is to transmit the pilot code S(z) prefiltered by -CL(z) from the RHC 
antenna and by CR(z) from the LHC antenna. After propagation through respective 
channels, the mobile terminal receives 

(-CL(z) -CR(z)+Cr (z) .CL(z) )S(z) = 0 

■if the prefiltering channel polynomials are indeed equal to the actual 

propagation channels. The pilot code would be received strongest at a virtual 

terminal location to which the propagation channels were -CL ( z ) # 

for RHC and CR ( z ) # for LHC . The propagation matrix for the intended 

signal and the pilot signal to the real mobile terminal and the virtual terminal 

location is thus 



CR CL 



-CL # CR # 



and is now square. The prefilter for transmitting X(z) to the real mobile 
and S(z) so as not interfere with the real mobile is the adjoint of this, i.e. 



CR # CL # 
-CL CR 

or in terms of the notation Pl(z) and P2(z), is 
Pl(z) P2(z) 
-P2 ( z) # Pi (z) # 

and the loopback signal received from the real terminal is 

(Pi ( z ) CR ( z ) +P2 { z ) . CL ( z ) ) X ( z ) + (Pl # (z)CL(z)-P2 # (z}CR(z) )S(z) 

Correlation with the independent symbol streams X(z) and S(z) after compensation 
of the loopback signal for the uplink channel yields 

Pl(z)CR(z) + P2(z)CL(z) = result of correlation with X(z) 
-P2 # (z)CR(z) + Pl # {z)CL(z) = result of correlation with S(z) 

These are now sufficient equations to determine CR(z) and CL ( z ) . 

The informational deficiency problem for the case of more antennas than 
terminals was solved above by imagining virtual or dummy terminals to exist 
for which the propagation coefficients are in some sense orthogonal to those of 
the real terminals, and to transmit pilot symbols streams to the dummy terminals 
that are statistically indepedent and preferably orthogonal. The additional 
row used to square up the 2x1 matrix to a 2x2 matrix in the above example may 
be seen to be the conjugate time-reverse of its own cof actor. 
The determinant of the 2x2 matrix is thus of the form 

CR.CR # +CL.CT.r 

which is the same as [C.C^]" 1 where C is the original 

2x1 matrix. The addition of the dummy terminals and propagation channels has not 
therefore affected the poles of the denominator of the overdimensioned problem. 



The gain of the overdimensioned solution therefore remains intact. 

The solution for three-base, two-mobile case may now be deduced as follows: 

Let the 3x2 propagation matrix be 

Cll(z) C12(z) C13(z) 
C21(z) C22(z) C23(z) 

The "dummy" terminal should have dummy channels forming a third row of 
polynomials that is orthogonal to both of the above. A 3 -element vector V 
perpendicular to two given 3 -element vectors corresponding to rowl and row2 
can be formed from the cross product of rowl and row2 , namely: 

V = rowl X row2 = C12 . C2 3 -C22 . C13 ; C13.C21~C11.C23; C11.C22-C12.C21 

The virtual channel should be chosen to the time-reverse conjugate of this, 
namely 

row3 = V # 

or C13 = C12 # .C23 # -C22 # C13 # 
C23 = C13 # .C21 # ~C11 # C23 # 
C33 - C11 # .C22 # -C12 # C21 # 

which may be verified also to be its own cof actor. Thus, as in the expansion of 
the 2x1 case to case, the expansion from a 3x2 overdimensioned problem to a 
3x3 completely determined problem using the above virtual channel construction 
does not alter the poles of the determinant of the problem and therefore 
preserves the efficiency advantages of the overde termined solution. 

In the. case of three base stations collaborating to form a coherent 
macrodiversity transmission to a single mobile, two non-polluting pilot-code 
transmissions must be added. When two vectors are given in a 3D space, a third 
direction perpendicular to them is unique; however, when only one vector in a 3D 
space is given, constructing another two to form therewith an orthogonal base is 
not not unique but invloves an arbitrary rotation about the given vector. 

For example , given channels 

Cll C12 C13 

a first dummy channel set 0 -C13 # C12 # 

for transmitting a first pilot code can be determined. Then the third row for 
transmitting a second pilot code can be determined to be equal to the 
time-reversed conjugate of its own cofactor, namely 



C31 - C12.C12 # + C13 # .C13 
C32 - -C11 # .C12 
C33 = -C11 # .C13 

However, other choices for row2 could be found which would then' result in other 
choices for row3 . All choices however preserve the denominator poles and thus 
the efficiency of the 3x1 overdimens ioning . 

In general a mobile terminal can be located within a triangle formed from the 
three nearest base stations, _ and therefore it is not necessary to extend the 
discussion to matrices of polynomials of dimension greater than 3. However, 



it may be deduced by extrapolation that the virtual Nth row in an Nx(N-l) 
overdimensioned solution would turn out to be the time-reversed conjugate of 
its own cof actor in the NxN matrix. Likewise in an Nx(N-2) problem, the 
(N-l)th row can be determined by setting any polynomial element to zero and 
making the remaining N-l polynomials the time-reverse-conjugate of "there own 
cof actors in the (N-l )x (N-l) matrix formed from the remaining rows and columns 
not comtaining the zero polynomial. Then the Nth row would be filled in as 
before . 

It has been described above how pilot codes necessarily transmitted by base 
stations to allow determination of the propagation channels from the base 
stations to a reduced number of mobile terminals can be prefiltered for 
transmission such that they cause no additional interference or pilot pollution 
to those mobile receivers, thus improving performance compared to the prior art. 



4. CLAIMS 



4. A method of radiating signals from a first number of transmitters helpful 
in the determination of propagation channel characteristics to a second 
number of receivers less than said first number/ comprising: - 

-transmitting wanted information signals to said receivers using said 
transmitters jointly such that the information signals are communicated 
most efficiently; 

-transmitting a third number equal to the difference of said first and 
second numbers of known pilot signals using said transmitters jointly 
such that said pilot signals cause substantially no interference to the 
reception of said wanted information signals. 
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